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Abstract 

We discuss and prove validity of the Hubbard- Stratonovich (HS) iden- 
tities over hyperbolic domains which are used frequently in the studies on 
disordered systems and random matrices. We also introduce a counterpart 
of the HS identity arising in disordered systems with "crural" symmetry. 
Apart from this we outline a way of deriving the nonlinear cr-model from 
the gauge-invariant Wegner k— orbital model avoiding the use of the HS 
transformations. 



1 Introduction 

A considerable progress achieved over the past two decades in understand- 
ing statistical properties of a single electron motion in disordered and chaotic 
systemsfl, 2] is mainly based on the nonlinear cr-model description. This con- 
cept was originally proposed in the context of disordered systems by Wegner [3], 
and further clarified and developed to a working tool in the papers by Schafer 
and Wegner [4] and Pruisken and Schafer [5]. The authors reduced the prob- 
lem of evaluating the disorder-averaged correlation functions of the resolvents 
(Green's function) of a Hamiltonian containing disorder to an effective deter- 
ministic field theory with a very peculiar underlying non-compact symmetry, 
which they called "hyperbolic". Original derivations [4] used a specific micro- 
scopic model, the Wegner's gauge invariant k— orbital model 1 . It was however 
soon understood, that the reduction is valid under much more general condi- 
tions. Loosely speaking, the nonlinear cr-model adequately describes physics of 
the one-electron motion for scales longer then the so-called mean free path. In 
the gauge invariant model the mean free path is effectively zero, which makes 
the derivation particularly transparent. 

In all the mentioned work the reduction of the microscopic disordered mod- 
els to the nonlinear cr— model made use of the so-called replica limit. The last 
essential ingredient - the idea of supersymmetry, i.e. the use of both commuting 
and anticommuting variables to avoid problematic replica method - was intro- 
duced in the theory by Efetov [1]. The resulting theory was successfully tested 
in an important non-perturbative limit of "zero" spatial dimension, where its 
predictions were shown to be identical to those following from the theory of 

1 Readers not familiar with the k— orbital model and its relation to nonlinear ct— model may 
consult a recent preprint [6], or Appendix A of the present paper for a short introduction. 
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large random matrices. After that within few years the supersymmetric version 
of the nonlinear er— model was accepted as a standard tool in condensed matter 
physics, and also proved to be very useful in other fields ranging from the theory 
of chaotic scattering [7, 8] to Quantum Chromodynamics[9]. 

Successful as it was in applications, the nonlinear a— model is derived from 
underlying disordered Hamiltonian along a mathematically subtle procedure. 
One of the standard non-trivial ingredients of the derivation is the exploitation 
of the so-called Hubbard-Stratonovich identity 2 : 

C N e-^ A ' 2 = JvR e -^-*TWM (1) 

where both R, A are square n x n matrices with entries Rij and Aij , respec- 
tively, and C/v is some A— independent constant. When we take these matri- 
ces to be complex Hermitian (or real symmetric), use the corresponding "flat 
differentials" VR oc Il"=i dRu Oj<j dlm[Rij]d~Re[Rij] as volume elements for 

Hermitian case, and DR oc H"=i dRu Yii<j dRij for real symmetric case, and 
integrate in the right-hand side over all degrees of freedom independently from 
minus to plus infinity, the integral amounts to a product of n 2 (resp. n(n+l)/2) 
one-fold Gaussian integrals, and the identity is completely trivial. The tricky 
point is that in the problem under consideration such a simple choice of the inte- 
gration manifold of matrices R is prohibited by some extra requirements which 
impose restriction on the structure of the matrices A and R. In particular, for 
the whole theory to be well-defined one has to ensure second (linear in A) term 
in the exponent in the right-hand side to be purely imaginary, by the very way 
one makes use of the HS identities. With a non- Hermitian choice of the matrices 
A,R it is already a rather nontrivial task. For an informal discussion of this 
point in a simple example see [10]; rigorous and comprehensive treatment can 
be found in the recent review paper by Zirnbauer[ll]. 

The difficulty was first encountered and successfully solved by Schafcr and 
Wegner. For the model which stems on the microscopic level from disordered 
Hamiltonians with broken time-reversal invarance (e.g. due to presence of a 
magnetic field) they suggested the following choice of the integration domain: 

R = XT f t + iP (2) 

where T G U(ni, n 2 ) is in the pseudounitary group of complex n x n matrices, 
n = ri\ + ri2, with integer n\ > 1, r&2 > 1. The inverse for such matrices is given 
by T _1 = LT^L, with L = diag(l ni , — 1„ 2 ) where 1„ stands for the identity 
matrix of size n, and stands for the Hermitian conjugate of T. The matrices 
P are Hermitian block-diagonal: P = diag(P ni , P n2 ) = an d the parameter 
A > is arbitrary. For the case of Hamiltonians respecting the time-reversal 

2 The original paper [4] does not use the name Hubbard-Stratonovich, but such a terminol- 
ogy became standard after Efetov's work[l] . 
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invariance the structure of the integration manifold is very much the same, with 
the pseudoorthogonal group 0(n\,n2) replacing the pseudounitray one. 

It is the manifold of matrices T that encapsulates non-compact ("hyper- 
bolic") symmetry of the problem, and at the next stage gives rise to the in- 
teracting Goldstone modes described by the nonlinear o— model. With such a 
nontrivial choice of the integration manifold the indcntity (1) loses its trans- 
parency, and its verification is a separate, nontrivial task. The idea of the proof 
[4] is to interprete the corresponding integral as going over a high-dimensional 
contour deformation of a simple Hermitian matrix. Accurate implementation of 
this argument can be found in the mentioned review paper by Zirnbauer [11]. 

Although the Schaier-Wegner parametrisation (2) of the integration manifold 
is admissible it was never much in use, virtually abandoned in favour of an 
alternative one due to Pruisken and Schafer[5]: 

R = f- 1 Pf, VR = d / j lH (T)dP 1 dP 2 A 2 [P} (3) 

where P = diag(P ni , P n2 ), with P ni and P n . 2 being real diagonal 3 , and T is 
again as in Eq.(2). The notation d/j,u(T) is used here for the invariant (Haar's) 
measure on the pseudounitary group, and for any diagonal n x n matrix B = 
diag(6i, 6 2 , . . . , b n ) we use the notation A[B] = Ili<j(^ — fy) for the associated 
Vandermonde determinant. The matrices A relevant in present context are 
always of the form A = A + L, where A + > is Hermitian positive semidefinite, 
and L is the "signature" matrix entering the definition of the pseudounitary 
group U(ni,n 2 ), see above. The Pruisken-Schafer type of parametrization for 
the integration manifold is used in the majority of applications (in particular, in 
the supersymmetric version of the theory, implicitly in [1] and explicitly in [7] and 
subsequent papers). Nevertheless, it seems that validity of the corresponding 
Hubbard-Stratonovich transformation Eq.(l) was never properly checked. In 
fact, for many years it was taken for granted[10] that a kind of "shift of the 
contour" argument should be valid for such a parametrization as well. The 
belief was however challenged by Zirnbauer (see criticism in[14], and recently in 
[11]) who revealed via a thorough analysis serious difficulties in implementing 
this type of argument for Pruisken-Schafer parametrization. This observation 
makes the situation uncertain and calls for further investigations. 

It is worth mentioning, that recently an alternative method of treating the 
hyperbolic symmetry was introduced which avoids the use of any variant of 
the Hubbard-Stratonovich integral[12], see a very informative discussion in [11]. 
This method can be successfully applied to gauge-invariant fc— orbital model, 
and we outline the corresponding derivation of the (bosonic) nonlinear er— model 
in the Appendix A of the present paper (see also [6]). Nevertheless, the Hubbard- 

3 The convergence of the integral (5) in the text below requires, in fact, infinitcscmal imag- 
inary parts to be included in variables P as P = diag(P ni — iO+l n , P„ 2 + «0+l n ). This shift 
is henceforth assumed implicitly. Similar remark is applicable also to the formula (25) in the 
"chiral" case. 
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Stratonovich transformations remains an important (frequently, the only avail- 
able) tool for majority of microscopic models, in particular for the popular 
model of band matrices, see [13] and discussions in [11], as well as for the case 
of random potential ("diagonal") disorder. 

The main goal of this paper is to provide a way of independent verification of 
the validity of the Hubbard-Stratonovich identity for the Pruiskcn-Schafer type 
of paramctrizations. We succeed in solving the problem in full generality for 
the pseudounitary case. As to the pseudoorthogonal case, the situation turns 
out to be more complicated and is in fact quite interesting. By considering 
explicitly the simplest nontrivial case n\ = 712 = 1 we will show that to save the 
identity one has to discard the properties for the factors entering the elementary 
integration volume to be positive and for the constant Cn to be real. Finally, 
we will discuss and prove the counterpart of the Hubbard-Stratonovich identity 
which arises naturally in the studies on disordered systems with special ("chi- 
ral") symmetry, and is also useful when investigating non-Hermitian random 
Hamiltonians. Some technical details are provided in the Appendices B and C. 



2 Hubbard-Stratonovich identities over Pruisken- 
Schafer domains. 

2.1 Pseudounitary case. 

We start with proving (1) for the general case of the pseudounitary Pruisken- 
Schafer domain R = T~ 1 PT parametrized by matrices T G [/(rii,n 2 ) and 
P = diag(Pi, P2) as in (3) . It is easy to show (see e.g. [12], Appendix A) that for 
strictly positive definite matrices A + > the corresponding matrices A = A + L 
always can be parametrised as A = T^ATJ 1 , where A = diag(A m , A„ 2 ) is 
diagonal such that A m > 0,A„ 2 < 0, and Ta is a pseudounitary matrix from 
U(ni, n 2 ). Starting with this form of the matrix A, we see that the right-hand 
side of (1) can be written as 

I { jPf{A) = J dPA 2 [P}e-^ Tlp2 J dn H (f) cxp{-iTrf- 1 Pf f A kf ^} (4) 

Now it is obvious that the integral is independent of the matrix Ta- Indeed, 
that matrix can be absorbed into T by the change of variables TTa — * T, 
exploiting the cyclic invariance under the trace and the invariance property of 
the Haar's measure: d/i^(T) = cIhh(TTa)- Second observation is that due to 
(block-)diagonal structure of the matrices P, A the combination in the exponent 



Tr 



T^PTk 



stays invariant when T is multiplied with an arbitrary unitary 

block-diagonal matrix of the form diag(Vi, V2), with V\ <G U(n\) and € U{n2). 
Thus, the combination in question does not change if we replace T with matrices 
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To € u(n^®u(n2) taken from the cosct space obtained by factorizing the original 
pseudounitary group U{n\, n 2 ) by its maximal compact subgroup U {ni)®U (n 2 ). 
As the result, the integral - up to a multiplicative constant - can be replaced by 
one going over the coset space rather than the whole pseudounitary group. This 
is a very pleasing fact, since such an integral has been recently shown[15] to be 
exactly calculable with help of the so-called Duistermaat-Hcckman localization 
theorem, generalizing a similar formula known for the unitary group[19]. The 
result of such integration is given, again up to a multiplicative constant, by 4 

J A[P]A[A] V ' 

Substituting (5) back to (4), we bring the latter to the form 

K AjAj / dPldh A[P] e " T " 2 det ^ iPuXl ^j=i det [ e " iP2iA2 C-=, 1+1 

(6) 

Now we observe that A[P] = A[Pi]A[P 2 ] U7=i U]= ni +i(Pu ~ and usc thc 
invariance of the integrand with respect to any permutation of the indices of 
integration variables in the set pu, i = 1, . . . , ni, as well as in the set p2j, j = 
ni + 1, . . . , n. Such an invariance allows one to replace each determinantal factor 
in the integrand with only one (" diagonal" ) contribution, multiplying the whole 
integral with the factor n\\n2\. Disregarding the multiplicative factors, we write 
the resulting integral as 

(i) o, _J_y iAdftAlPle-i^e-'^.^^e"'^^ (7) 

and further use the following well-known identity 5 

' dPA[P] e -^Tr^±iTr[PA] a A[ A] e -|TrA= 



valid for any real diagonal matrices P, A, . In view of TrA 2 = Tr^4 2 the lat- 
ter result shows that the original integral (4) is equal, up to a multiplicative 
^.—independent constant to e~^ TA . We thus proved the required Hubbard- 
Stratonovich identity for matrices A = A+L, with A + > 0. 

What remains to be shown is how to incorporate the case of positive semidef- 
inite A + > in the above scheme, which is crucial for applications. The problem 
is that if the matrix A + has zero eigenvalues, the corresponding matrix A = A + L 
may be not T— diagonalizable. 6 To this end introduce a parameter e > and 

4 See the previous footnote on convergence of this integral. 

5 This formula can be for example derived starting from (1) for Hermitian matrices R, A, 
diagonalising R by unitary transformation, and performing the unitary group integration 
by the Itzykson-Zuber-HarishChandra formula[19]. See Appendix B for a proof of a similar 
expression in "crural" case. 

6 The simplest relevant example is A = ^ j 1 ) ' w ' 1 ' c ' 1 sc l uares to zcro: A 2 = 0- 
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consider the integral: 

y*X)7^ c -2Tr(R+i e L) 2 -iTrAR . = e i e 2 Trl„ y c -lTr R 2 -iTr (A+eL)R _ ^ 

From A = A + L, where A + > immediately follows that A f = A + eL = 
(A + + el n )L = A f __ + L, where A e _ + is already positive definite: A t _ + = A + + 
el n > 0. But such matrices A e are always T-diagonalizablc, and the above- 
given proof of the Hubbard-Stratonovich identity retains its validity. Therefore, 
the integral is calculated as 

e ie 2 Trl n y p^ e -±TrR 2 -iTr(A e R) = g i e 2 Trl ne - iTr (A e ) 2 _ g— eTr LAg— iTr A 2 

resulting in the e— modified version of the Hubbard-Stratonovich identity: 

J Pft e -5 T r(R+i£L) 2 -iTrAR = e _ e Tr LAg-±Tr A 2 ^ 

For e > this identity holds uniformly in A, including the case of non-diagonalizable 
matrices A. 



2.2 Pseudoorthogonal case: ri\ = = 1. 

The difficulty of proving (1) in the important case of general pseudoorthogonal 
group 0{n\, rz 2 ) is due to lack of integration formulas similar to Eq.(5) for cosets 
\ ■ Under these circumstances we will restrict ourselves by the first 
non-trivial case n\ = n 2 = 1, which proves to be already very informative. The 
matrix P in this case is 2 x 2 diagonal: P = diag(pi,p 2 ), and the matrices Tq 
can be explicitly parametrized in terms of the variable 6 e (—00,00) as Tq = 

( srnh0 cosh# ^ en tnc 2x2 matrices R = Tq^PTq we are integrating 
over in Eq.(l) are explicitly given by: 

/ Ei±P2 + £±_P2 cosh2 e £i_£2. sin h2# \ 

^ _£i_P2 sinh2e P1+P2 _ pi_£2 cosh2 6» J (W) 

As we already mentioned, the matrices A must be of the form A = A + L, 
where L = diag(l,— 1). We restrict ourselves in this section only with real 

symmetric matrices A + > for simplicity, i.e. A + = j fll a ] > 0. Hence 



a a 2 



a — a 2 



A = ( ai a I , with the constraints 



ai > 0, a 2 > 0, \a\ < ^/aia 2 (11) 



6 



Naively, one may expect the choice of the volume element on such a manifold 
in the form dR = |pi — P2\dpi dp2 d6 to be "natural" . We will however see below 
that taking such a choice we end up in a trouble: the corresponding Hubbard- 
Stratonovich formula (1) does not hold its validity any longer. Instead, the 
correct choice of the "volume element" in our case turns out to be: 

dR= (p 1 -p 2 )d Pl dp 2 d8 (12) 

This expression is not sign-definite any longer, and changes sign at the line 
Pi = P2] we will comment on this point shortly later on. 

After having specified all ingredients of the right-hand side in (1) we can 
write down the corresponding integral explicitly as 

/OO poo 
dpi dp 2 (pi-p 2 )e-^ +pl) - i ^ pi+P2){ai - a2) (13) 
-OO J —OO 

dQg~i a cos h 29—ij3 sinh 20 



where we introduced the shorthand notations a = \{a\-\-a2){pi — P2), = a(pi — 
p 2 ). We note that in view of (11) holds 0/a = ^ < 7 = < 1. Therefore, 

we can parametrize (3 = usinh-0, a = ucosh-!/;, where u,ip are real parameters. 
Then the combination entering second exponent in (13) can be rewritten as 
a cosh 29 + (3 sinh 29 = w cosh (2$ + tp). Finally, introducing \i = 29 + ip as 
integration variable, we see that the integral over 9 can be explicitly calculated 
as [16] (a): 

-J dve- mcosh ^--{Y Q (\u\)+2 S gn(u)M\u\)} = K (iu) (14) 

where sgn(u) = ±1 depending on the sign of the variable u, and Jo (z), Y (z), K (z) 
are Bessel, Neumann and Macdonald functions of zero order, respectively [16]. 
In our case: 

U = (pi-P2)s a , s a = \ ( ai + a2 ) -a 2 (15) 



Substituting the result of integration back to (13) and changing to integration 
variables p± = (pi±p 2 ), one can easily perform the Gaussian integral over p + 
and obtain 

/OO 
dp-p-e-^- [Y (\p_\s a )+i S gn(p_)J Q (\p_\s a )] (16) 
-oo 

This is exactly the point where we can most clearly see the necessity of omitting 
the absolute value sign in the measure (12). Indeed, had we maintained the 
modulus in the above integral, the second term in the integrand, being 



7 



odd in p-, would vanish and the remaining integral would be that containing 
the Neumann function. Although it is well-defined, it could not produce the 
structure necessary for the validity of the identity Eq.(l). In contrast, when 
the factor p_ in the measure does not contain the absolute value, it is the first 
term which vanishes, and the second term can be straightforwardly integrated 
by using the identity[16](b): 



exactly as required by the Hubbard-Stratonovich identity Eq.(l). 

The considered example (and also one in the next section) makes it clear that 
the only consistent way of ensuring the validity of the HS transformation over 
hyperbolic domains is to require absence of the absolute value of non-positively 
defined factors in the elementary volume. From geometric point of view the in- 
tegration domain in hyperbolic case consists of several (in the simplest case two) 
disconnected pieces. One may notice that the factor in the elementary volume 
(12) changes sign precisely when passing from one such piece to a different one, 
being sign-constant within each piece. Persistence of this structure for n > 1, as 
well as finding a comprehensive proof valid for any n in the pseudoorthogonal 
case remains an interesting open question deserving further attention. 

2.3 Chiral variant of the Hubbard-Stratonovich identity 

Last decade new symmetry classes of random Hamiltonians attracted a lot of 
interest due to numerous applications in various branches of physics, see [17] for 
discussion and basic references. In particular, the class of Hamiltonians with chi- 
ral symmetry is pertinent for analyzing properties of Dirac fcrmions in random 
gauge field background, and found applications in Quantum Chromodynamics 
[9], as well as in condensed matter theory, see e.g [18] for more references and 
discussion. When reducing analysis of such systems to the relevant nonlinear 
u— model, one encounters the following variant of the Hubbard-Stratonovich 




yielding, with b = 1/4, c = s, 



•a 



j(po) _ _ 4i7] .3/2 e -i(a 1 -a 2 ) 2 -I(a 1 +a 2 ) 2 +a 2 



(17) 



identity: 




In the simplest case the two involved matrices are related as A 1 = B e GL(n, C) 
, i.e A is an arbitrary complex matrix, and one can make a natural choice of the 
integration domain R[ = R 2 € GL(n,C), with elementwise "flat measure" on 
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it. The linear in A term in the exponent is then purely imaginary as required, 
and the identity 



r(ch) 
1 HS 



(A, it) = C N e- TT ^ = J VRVltf e-Tr&R-iTrl&A+^R] (19) 



follows from the integral in the right-hand side being the standard Gaussian 
one. 

In the applications, however, the case of two unrelated Hcrmitian positive 
scmidefinite matrices A' = A > 0, B'=B>0 emerges as well, and the extra 
convergency arguments necessitate to make the choice for the integration domain 
to be compatible with that property (see below). Again, ensuring the pure 
imaginary nature of the exponent and subsequent verification of the Hubbard- 
Stratonovich identity (18) for such a domain is a non-trivial task which was not, 
to the best of author's knowledge, yet accomplished in full generality. 

Our strategy in this case will be informed both by our experience with the 
pseudounitary and pscudoorthogonal cases. Again restricting ourselves in this 
section with A > 0, B > we observe that any pair of Hermitian, positive 
definite n x n matrices A, B can be parametrized as (see [18] appendix A) 



ft 

1 A 



aTX 1 



A = T A aT\, B 
where a — diag(a!, . . . , a„) > 0, T A € 



GL(n,C) 



U{1) x ... x u(iy 



(20) 
(21) 



that is the matrix T4 is a general complex with real positive diagonal entries. 
This suggest an idea to parametrize the integration manifold as 



Ri 



1 „ r „n -1 
T t p T 



R 2 = fPf f 



(22) 



in terms of a real diagonal matrix —00 < P = diag(pi, . . . ,p n ) < 00 and a 
general complex matrix T G GL(n,C). Guided by our previous experience, the 
volume element in the new coordinates is chosen by: 



VR\VR 2 cx Y[pid Pl J] {pi - pl) 2 dfi H (f, ft) 



(23) 



1=1 



l<m 



where d/x_y(T, Tt) is the invariant Haar's measure on the group GI(n,C). Note 
that despite the domain of integration being — 00 < ^ < 00, i = 1, . . . ,n the 
volume element contains factors pi rather than \pi\ as one may naively expect. 

Substituting such a parametrization into the integral (18) and using the 
cyclic invariance under the trace, we have Tr[_Ri_R 2 ] = TrP 2 and the T-dependcnt 
term in the exponent is given by 



Tr[i?ii + BR 2 ] = TrP 



f^f A af A f + f- 1 (f A ^j 1 aTV 1 (ft) 
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Now we change variables T\T — > T and exploiting the invariancc of the mea- 
sure d/z#(T, T 1 *) = dfiH(T,T^) satisfy ourselves that we need to deal with the 
following group integral: 



TeG\(n,C) 



dfj, H (T,Ti)exp-iTrP 



T'aT + T~ a ( T 



(24) 



Again, due to diagonal structure of the matrices P,a the integrand is not 



changed if we replace matrices T € GI(n,C) with Tq € 



GL(n,C) 



This 



(7(l)x...xt/(l) • 

fact provides us with the possibility of exploiting one more integration formula 
discovered in [18] 7 : 



/ 



Gl(n.C) 



d/z (T 0) T T )e 



- 4 Trp[f tafo+f- 1 a(f ( t)- 1 ] det [K {2i Pl a 



3>H,j = l 



[U(1)X.-.XU(1)] 



A[P 2 ]A[a 2 



(25) 

where d/ioiTo,^) = dT^dT^ det TbTj 2 is exactly the volume element on 

the coset space manifold. Substituting this result back to the integral (18), 
cancelling one of the Vandcrmonde factors coming from the volume element 
(23), and again exploiting the invariance properties of the integrand with respect 
to permutation of integration variables in the set p\ , . . . , p n (cf. discussion 
after (6)), and disregarding A— independent multiplicative factors, we bring 
the integral of (18) to the form 

4f(AB) cx — — / Y[dp lPl A[P 2 ]e-2:^ l[K (2 t p l a l )(26) 



i=i 



dpp e- p2 p 2si K (2ipai) 



(27) 



where we used the expansion of the Vandcrmonde determinant as A[P 2 ] = 
Ss Q (—1) niLi P 2si m terms of the sum over n\ permutations S a = (si, s 2 , . . . , s„) 
of the index set (1, 2, . . . , n), with (— l) Sa standing for the parity of the permuta- 
tion. Now we substitute the formula (14) for the Bessel function of the imaginary 
argument to the above expression, and see again that the term with Neumann 
function is multiplied with the odd in p factor and vanishes upon integration, 
yielding 



J<#(i,B)« ' 



A 2 [a 2 ] £ 



n r „ 

(-I)S-IJ (_ i7r) / 

1 Jo 



{-in) I dpp e-P p 2si J (2p ai )(28) 
o 



7 Although the formula itself is correct, its derivation in Appendix B of [18] was not accurate 
enough. For this reason we include the outline of the correct derivation in the Appendix C to 
this paper. 
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where we used sign(ai) = 1. Finally we invert the operation of the Vandermonde 
determinant expansion and show in the Appendix B that the resulting integral 
can be evaluated as : 

exactly as required by the Hubbard-Stratonovich identity (18). 
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Appendix A: From /c-orbital model to nonlinear 
o— model without Hubbard-Stratonovich identity. 

Consider N x N random Hcrmitian matrix H, with all its entries Hi m for 
I < m being independent random Gaussian with the variance (Hi m Hi m ) = 
Jim > 0. Let N — rk, with r and fc being integers. Subdivide the index 
set / = 1,2, ...,7V into r subsets Ii,I 2 ,...,I r each having exactly k elements: 
1% = {{i — + 1) (i — + 2, (i — l)k + k} and consider the variances J; m 
such that 

{J/k , if(Z,me/») 
V/k 2 , if(Z e Ii,m£ or (I e Ii-um G h) (30) 

otherwise 

This defines the Wegner's gauge-invariant fc— orbital model[4] of r "sites" ar- 
ranged in the one-dimensional lattice d = 1. Each "site" is a block Ii represent- 
ing a group of k orbitals, with J/k standing for intragroup coupling and V/k 2 
representing the coupling of neighbouring groups. The matrix H can be visu- 
alized as having a (block-)banded structure made of fc x fc blocks, with nonzero 
entries concentrated in a band of the widths oc fc around the main diagonal. The 
couplings are scaled in a way ensuring correct behaviour in the limit of infinitely 
many orbital fc — > oo, which is assumed at a later stage (see also closely related 
model in[23]). Generalization to lattices of higher spatial dimensions is obvious. 
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Our goal here is to demonstrate that the method developed in [12] for random 
matrix problems without underlying lattice structure can be straightforwardly 
applied to arrive to a nonlinear a— model from the H defined above. We start 
with considering the simplest nontrivial object, the negative integer moment of 
the modulus of the spectral determinant of H: 

Z n , n {E, n) = (l dot (E + iv/k - H)\- 2n ) H (31) 

where E and 77 > are real and brackets stand for the ensemble averaging. Our 
first goal is to show that it has the following integral representation: 

Z n , n (E,r)) = Const x J d qi ... J dq r e~ k ^i=i £(<5l) (32) 

fildct^] 2 " 
where L = diag(l„, -1„) 

C(q) = ^Trf - iETrq - Tr In q (33) 

and the integration goes over the positive definite Hermitian 2n x 2n matri- 
ces qi > 0. This representation is exact (no approximations) with the only 
restriction k > 2n. 

Here wc outline the main steps employed to derive Eq.(32). For every 
I = 1, TV introduce two n— component vectors $; !+ , with complex COin- 
^l) c (n) nnA n i^ n , 10 „ ^ n 1T1. — I '>+ 

Representing the inverse determinants as Gaussian integrals over $;,± we have 



poncnts Si ±, ...,S , ; ^_ and also use the notation = f J and rjk = rj/k. 



r- r- N 

Z ntn (E, V )oc J d*!...y d9 N exp{Y^*\ (iEL -n k l 2n ) *!> (e^E.™"'-"*^*™^ (34) 

Now for any I = 1, . . . , N we introduce 2nx2n Hermitian matrices Qi = ^i®^], 
so that the result of the (gaussian) averaging for any variances Ji m can be written 

as 



I H 

For the particular choice of the variances Eq. (30) we can further write 
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where we denoted Q s = Qk(i-i)+ s - Introduce now the set of matrices cji = 
J2s=i f° r * = 1' •••' r anc ^ rewr itc the integral as 

Z n , n (E, V )ot J ^ 1 ...^ r X(g 1 )...X(q r )e-* Tr SL 1 (^ i ) 2 -i&E[: i lTr (^ £ )(^^) (35) 
where 

T(q) = J d*i...d* fc (S ^-2*.®*^ ex p{XJ*s (*^-»?fcl2n) *.|(36) 

Representing the (5— function factor as the Fourier transformation over the Hcr- 
mitian 2n x 2n matrix K: 

and taking into account Tr if <g) = <3>tX\|/ s W e find that the 'J— integrals 

are gaussian (and convergent in view of r/k > 0) and so when performed explic- 
itly yield the factor: 



det 



K - EL- ir) k l 2n 



i -k 



We immediately see that the resulting integral over K is precisely one of Ingham- 
Siegel type calculated in [12] 8 . For k > 2n it gives: 



T{q) = 0(q) (detq) 



^k-2n iTrq(EL+i Vk l 2 „ 



(37) 



where the factor 9 (q) is non-zero for positive definite matrices and zero oth- 
erwise. Substitute this expression back to Eq.(35) and finally change q — > k q 
This immediately produces the formula eq.(32). 

Now change the integration variables in Eq.(32) from q^ for i = 1, . . . ,r to 
the matrices &i = qiL parametrized as &i — T^ 1 P i T i , with Tj e U{n,n)/U(l)® 

... ® (7(1) and Pi = diag(pg, ...,pS), such that f>$ > O.pg 5 < 0. 

Correspondingly for each i — 1, r the integration measure is given by d&i cx 
A 2 (Pi)dfi(Ti)dPisdP2,i and the above integral assumes the form: 



Z ntn (E,rj) = Const x J da x ... J d& r e k ^ = i c(p ^ 

r 1 

J e --rZ)[Ji TrS '''* i + 1 " ? 'Z][=i Tr * i ^ 



(38) 



n 



dot P, 



8 A mathematically-minded reader who may dislike a little bit frivolous manipulations with 
S— functions, may wish to perform the derivation by exploiting the "integration theorem" 
proved in [15]. See also an alternative way of derivation in [6]. 
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where we denoted 

C(P) = ^-TrP 2 - iETrP - Tr In P (39) 

Now it is evident that in the limit k — > oo all the matrices Pi will be fixed 
by the (unique) saddle-point value Po minimizing the "action" £(P) and given 

by 

At the next step one should accurately perform the integration of Gaussian fluc- 
tuations around the saddle-point ( presence of the Vandcrmonde determinants 
makes the pre-exponential factors vanishing at the saddle-point for any n > 2). 
The resulting expression will clearly have the following structure: 

Z n , n (E, V ) oc J dQ^.J dQ r {...} e-^EI^^S^-^EL^L (41) 

where the integration manifold is parametrized by the matrices Qi = T^ 1 LT i . 
Here we used dots to denote pre-exponential factors which may possibly arise 
and also set the energy parameter E = for simplicity. This is exactly the 
lattice version of the nonlinear a— model introduced by Wegner and Schafer[4]. 
The index i = 1, . . . , r numbers the lattice sites and for r — > oo in the lattices 
of dimensions d > 2 the effective coupling constant controls the transition 
from localised to extended states in the underlying Hamiltonian H . 

The saddle-point procedure at k — > oo and subsequent manipulations are still 
to be done rigorous by strict mathematical standards, see recent progress and 
discussions of related issues in [6, 24] . Another interesting problem is how to 
include anticommuting degrees of freedom in the above derivation. Technically 
this can be done following various methods, and will be discussed elsewhere [25], 
see also [11]. 

Appendix B: proof of the formula (29). 

The starting point of the proof is the identity (19). Introduce the singular value 
decompositions of the general complex matrices R, A as: 

R = UPV\ P = diag(pi, . . . ,p n ) > 0, A = U A aV\, a = diag(oi, . . . ,a n ) > 

„ . . . ( 42 ) 
where U, Ua, V, Va are n x n unitary, with the associated Haar's measures 

d[i(U) and dfi(V), respectively. The volume element VRVw in new coordi- 
nates associated with the singular value decomposition is given by VKDW oc 

nr=i pidrn n l<m (Pi - pI^mu) <wn 
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Introducing U = U^Ua and V"t = V\V , and using the invariancc of the 
Haar's measures: dfi(U) — d/j,(U), dfi(V) = d/i(V), we rewrite the integral in 
the right-hand side of (19) as: 

J^fl dpan A 2 [P 2 } e" J d»{U) J d p<y )e -n±[H™W+*m\ (43) 

The integral over two unitary matrices in this expression is well-known from the 
papers by Guhr and Wettig [20] and Jackson et al.[21]: 

<^\Pli • • • iPn) L ^\ a li ' ' ' i u nl 

Now substitute (44) into (43), and after cancelling one Vandcrmonde factor 
observe that due to the invariance of the integrand all n\ terms of the deter- 
minant made of Bessel functions yield identical contributions. Disregarding 
the multiplicative constants one can therefore replace that determinant with 
JI™ Jo(2piO,i), so that the right-hand side of (19) takes the form 

1 r°° n 2 ™ 

— - / \JdpipiA[P 2 }e-^ p ' \{jv{2p m ) (45) 
\ Jo i=1 l 

On the other hand, the left hand side is exp —TtA^ A = e~ S;=i a ' which proves 
the formula (29). 

Appendix C. Matrix Macdonald functions asso- 
ciated with integrals over complex matrices 

The integral (44) can be looked at as certain matrix Bessel function that cor- 
responds to Itzykson-Zubcr-like integrals over unitary matrices, see details in 
Guhr and Wettig, Guhr and Kohler [20, 22]. Below we consider the matrix 
Macdonald functions that are associated with integrals over arbitrary complex 
matrices: 

/ dM f,ft )e -^(^[^ ft +(- t )" 1 ^- 1 ]) = const Af d T [Ko{x :H ] !T 3 - n 2 , 

(46) 

where Xd and Yd are two positive definite diagonal matrices: 

X d = diag(zi,:r2, • • • ,x n ) > 0, Y d = diag (j/i, y 2 , ■ ■ ■ , y n ) > (47) 

The integration in (46) goes over T e GI(n,C)/U(l) x ... x U(l), i.e. complex 
matrices with real positive diagonal elements, and d/x(T, T^) is the correspond- 
ing measure. Note, however that the integrand is not changed if we consider 
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f e GI(n,C) rather than restricting it to the coset space as above. Therefore, 
if we use the invariant Haar's measure d/x_y(T) for the full group, the result of 
integration will be the same up to a constant factor. 

Let us look at the integral in (46) as a function &(Xd, Yd). We start with 
considering a pair ofnxn Hermitian positive definite matrices X^ X \X^, and 
another such pair A, B. Introduce the Laplace operator D X (i) x < 2 ) acting on 
such matrices as: 



d 2 



Then we construct a function W(X^\ X^ 2 \ A, B) with the property 

£>x(D,x(2) W(X {1 \X {2 \A, B) = Tr(AB) W{X^\ X^ 2 \ A, B) 
In particular, the following function 



W{X {1 \X {2 \A,B) =exp 



(48) 



(49) 



(50) 



satisfies the Eq.(49) as can be checked by direct calculations. Let us now use a 
possibility [18] to parametrize A = TyYdT Y and B = (T^-^Tf 1 in Eq.(49), 
as well as X^> = T x X d T x ,X^ = (T x )~ 1 Xd(T x )~ 1 . Here both T x , Ty belong 
to the coset space GI(n,C)/U(l) x ... x U(l). 
We obtain: 



D X (i) X (2) exp--Tr 



X^T Y Y d Tl + {Tl)- 1 Y d T~ 1 xW 



Tr(y d 2 )exp 



^Tr (X d [f Y d f t + (f^Ydf - 1 ] ) 



(51) 



where we introduced the matrix T = T x Ty belonging to the group GI(n,C). 

Now we use that (i) the integration over complex matrices Ty commutes 
with the Laplace operator D X (i) X (2) and (ii) for any fixed Tx the result of 
integrating the right hand side of (51) over Ty does not depend on T x due to 
the possibility to use the invariant measure d^niTy) = duniT). We therefore 
conclude that the matrix function Q(Xd,Yd) defined by the integral Eq.(46) 
satisfies the following differential equation: 



D xW>xW *{X d ,Y d ) = Tr(Y 2 MXd,Y d ) 



(52) 



To derive the explicit formula Eq.(46) for the matrix function <f>(X d ,Y d ) we 
apply the method proposed by Guhr and Wettig [20] duly modified. We no- 
tice that when passing from X^\X^ to the "angular" coordinates T x and 
"radial" coordinates X d the radial part of the Jacobian is given by J{X d ) = 
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A 2 (x\ 1 . . . , x\ ) Y\ x i- This is enough to ensure that the radial part Dx d of the 

i=l 

Laplace operator D X (i) X m must have the following expression: 

n n 

Dx^y^Y, 9 ^^ J{x)=A 2 {x\,...,x 2 n )J{x i (53) 
J[x > i= i i= i 

Guhr and Wettig noted [20] that the radial part Dx d of such form is, in fact, 
separable. It means that for an arbitrary function f(xi,X2, ■ ■ ■ , x n ) the following 
identity holds: 

f(xi, ...,x n ) = 1 / d 2 1 9 \ f{xi,...,x n ) 

Xd A(xl...,xl) A(xl...,xl) f^\dxl + x k dxj A(xl...,xl) 

(54) 

The separability of the operator Dx d enables us to solve the differential equation 
Eq.(52) and to prove formula (46). The remaining parts of the proof given in 
the appendix B of [18] hold its validity without any modifications. 
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